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0.  Abstract.  »  t 

T\  ■_  .f  r-'.r'  / 

>-W»  address  the  problem  of  finding  a  minirpdm  weight  base  B  of  a 
matroid  when,  in  addition,  each  element  of  the  matroid  is  colored  with  one  of 
m  colors  and  there  are  upper  and  lower  bound  restrictions  on  the  number  of 
elements  of  B  with  color  i,  tor  i  -  1 , 2, ....  m.  This  problem  is  a  special  case 
of  matroid  intersection.  We  present  an  algorithm  that  exploits  the  special 
structure.  When  applied  to  the  weighted/bipartite  matching  problem,  our 
algorithm  has  complexity  0(|E|  |V|  +  |V|*  log  |V|).  When  applied  to  the  problem 
of  finding  a  minimum  weight  spanning  tree  with  degree  restrictions  on  the 
nodes  of  a  stable  set,  it  has  complexity  0(|  V|?).  In  both  cases,  V  denotes  the 
node  set  of  the  underlying  graph,  and  E  denotes  its  edge  set.  ‘/We  also  discuss  ^  • 
a  new  relaxation  for  the  traveling  salesman  problem. 

!  <S 

I 

.  J  <•  'VV  N 

1.  Introduction. 


Let  E  be  a  finite  set,  M  a  matroid  of  rank  r  defined  on  E,  and  w:  E  ->  R  a 
weight  function.  The  weight  of  a  subset  F  c  E  is  defined  by 

w(F)  -  X{w(e) :  e  e  F}.  Let  E1 ,  E2 . Em  be  a  partition  of  E  into  nonempty 

subsets.  To  each  Ej  assign  two  integers  lj  and  Uj  (lj  £  Uj).  We  address  the 


for  i  - 1, 2, ...,  m. 


problem  of  finding  a  base  B  of  M  which  will 
minimize  w(B) 

subject  to 

lj  £  |B  n  Ej|  £  Uj 

This  problem  is  denoted  by  (P). 

It  is  known  that  the  set 

B*-{XCE:  lj£|XnEj|£Uj,  i-1,2,. 

is  the  family  of  bases  of  a  matroid  M',  which  has  rank  r 
called  a  generalized  partition  matroid.  Thus,  (P)  is  the  problem  of  finding  a 
minimum  weight  common  basis  of  M  and  M',  if  one  exists. 

This  problem  has  been  solved  in  polynomial  time  by  Edmonds  ([El]  and 
[E2]).  The  algorithm  presented  here  takes  advantage  of  the  special  structure 
of  the  matroid  M'  and  works  on  an  auxiliary  digraph  with  m  nodes,  instead  of 
the  n  -  |E|  nodes  required  by  the  general  matroid  intersection  algorithms. 


.,  m,  |X|  -  r) 

Such  a  matroid  is 
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This  leads  to  a  complexity  of  0(r(nc  +  m  log  m)),  where  c  is  the  complexity  of 
circuit  finding.  By  circuit  finding  we  mean  the  following:  Given  independent 
set  I  and  e  e  E  - 1,  list  the  elements  of  the  unique  circuit  in  I  +  e,  or  show 
that  I  +  e  is  independent.  This  complexity  is  to  be  compared  with  the 
0(nr(r  +  c  +  log  n))  required  by  the  general  algorithm  (see  Brezovec, 
Cornu6jols,  and  Glover  ([BCG])). 

To  make  this  paper  self-contained,  we  state  without  proofs  some  of 
the  results  of  [BCG].  This  is  done  in  Section  2.  The  algorithm  for  solving  (P) 
is  presented  in  Section  3. 

Finally,  in  Sections  4  and  5,  we  provide  two  examples  of  how  our 
algorithm  may  be  applied  to  classical  combinational  optimization  problems. 
For  the  weighted  bipartite  matching  problem  on  a  graph  G  -  (V,  E),  we  get  an 
algorithm  of  complexity  0(|V|  |E|  +  |V|2  log  |V|),  which  is  the  best  known  for 
weighted  bipartite  matching  (see  Fredman  and  Tarjan  [FT]).  Note  that  the 
general  matroid  intersection  algorithm  only  gives  a  bound  of  0(|E|  |V|2). 

In  section  5,  we  consider  the  case  where  M  is  the  graphic  matroid 
associated  with  a  graph  G  ■  (V,  E).  Assign  each  edge  of  G  a  weight,  and  let 
V* « {uj,  v2, ....  vm.j)  denote  a  stable  set  of  G  (a  set  of  pairwise  nonadjacent 

nodes).  Let  color  class  i  consist  of  the  edges  incident  with  node  \)j  for 

i  -  1 , 2, ....  m-1 ,  with  color  class  m  containing  all  remaining  edges.  A 
minimum  weight  spanning  tree  of  G  must  be  found  with  the  property  that,  for 
i  *  1 , 2, ....  m,  the  number  of  edges  with  color  i  belongs  to  the  interval  [lj,  uj]. 

This  is  of  complexity  0(|V|^).  Note  that,  here  as  well,  the  general  matroid 
intersection  algorithm  only  gives  a  bound  of  0(|E|  |V|2).  If  lj  -  uj  -  2  for 

i  -  1 , 2, ....  m-1 ,  the  complexity  can  be  further  reduced  to  0(m  |V|2),  and  such 
degree  restricted  spanning  trees  can  be  used  to  strengthen  the  1  -tree 
relaxation  of  Held  and  Karp  for  the  Traveling  Salesman  Problem.  We  also 
show  how  this  relaxation  can  be  further  strengthened  and  generalized  to  the 
case  where  V*  is  not  a  stable  set.  In  this  case  the  color  classes  are  defined 
in  terms  of  net  stars  instead  of  stars. 


2.  Preliminary  results 

Let  M  be  a  matroid  defined  on  the  element  set  E,  and  let  B  be  a  base  of 
M.  Consider  p  c  B  and  p'  c  E  •  B  such  that  |p|  -  |p’|.  We  say  that  (p,  p’)  is  a 
B-swap  if  B  -  p  +  p'  is  a  base  of  M.  We  call  m(p,  p')  a  matching  if  it 
represents  a  one-to-one  mapping  of  p  onto  p';  m(p,  P')  is  called  a 
B-matching  if  every  (e,  o’)  in  m(p,  P')  is  a  B-swap. 

Remark  1.  [BCG,  Lemma  2].  Let  B  and  B'  be  two  bases  of  M,  and  let  p  -  B  -  B' 
and  p’  -  B'  -  B.  Then  there  is  at  least  one  B-matching  of  (p,  p’). 

Remark  2.  [BCG,  Lemma  3].  Let  B  be  a  base  of  M,  and  consider  p  c  B  and 
P'CE-B  where  |p|  -  |p'|.  If  m(p,  P’)  is  a  B-matching  but  (p,  p')  is  not  a 
B-swap,  then  there  also  exists  a  B-matching  m’(p,  p')  different  from  m(p,  p'). 

Consider  two  matroids  M1  and  M2  defined  on  the  same  finite  element 
set  E,  and  let  I  c  E  be  a  k-intersection,  that  is,  let  I  be  independent  in  M1  and 
M2  such  that  |I|  «  k.  Let  w:  E  R  be  a  weight  function.  The  problem  of 
finding  a  minimum  weight  k-intersection  is  denoted  (Pk).  To  solve  (Pk), 

define  a  bipartite  auxiliary  digraph  G(I)  *  (V,  V,  A)  as  follows.  Construct  a 
node  in  V  for  each  eel  and  a  node  in  V  for  each  e'  e  E  - 1.  For  each  e  e  I,  we 
construct  an  arc  (e,  e’)  with  weight  w(e,  e’)  -  w(e')  -  w(e)  provided  I  -  e  +  e' 
is  independent  in  M1 .  Similarly,  for  each  e'  e  E  - 1,  we  construct  an  arc  (e\  e) 

with  weight  w(e\  e)  ■  0  provided  I  -  e  +  e'  is  independent  in  M2. 

Let  C  be  a  dicycle  in  G(I).  We  define  the  weight  of  this  dicycle  to  be 
the  sum  of  the  weights  of  the  arcs  in  the  dicycle.  Thus,  w(C)  -  w(L')  -  w(z), 
where  I  and  r  represent  the  nodes  of  C  in  V  and  V\  respectively. 

Theorem  1.  [BCG,  Theorem  3].  Assume  I  is  a  k-intersection. 

(i)  I  is  optimal  for  (Pk)  if  and  only  if  there  are  no  negative  weight 

dicycles  in  G(I). 

(ii)  Let  C  be  a  negative  weight  dicycle  in  G(i)  with  no  negative  dicycle 
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on  a  subset  of  its  nodes,  and  let  £  and  £'  be  the  nodes  of  C  in  V  and  V', 
respectively.  Then  r  - 1  -  £  +  r  is  a  k-intersection  such  that  w(r)  <  w(i). 

This  theorem  leads  to  an  algorithm  for  solving  (Pk),  given  an  initial 

k-intersection  I:  use  G(I)  to  find  an  improved  solution;  continue  until  a  G(i) 
is  found  which  has  no  negative  dicycles. 

We  modify  this  technique  to  develop  a  dual  algorithm.  Consider  a  set  A 
of  k  artificial  elements,  where  EnA-0.  Let  F  -  E  +  A  and  define  two 
matroids  M1  (F)  and  M2(F)  on  the  element  set  F  as  follows.  For  i  -  1 , 2, 1  c  E 

and  J  c  A,  the  set  I  +  J  is  independent  in  Mj(F)  if  and  only  if  I  is  independent 

in  Mj.  The  problem  (Pk)  relative  to  these  new  matroids  is  denoted  (Pk|  F).  In 

order  to  define  (Pk|  F)  completely  we  need  to  assign  weights  to  the  artificial 

elements.  By  giving  them  large  negative  weights,  we  can  guarantee  that  I  «  A 
is  optimal  for  (Pk|  F). 

Define  the  digraph  Sp(l)  from  G(l)  by  splitting  one  of  the  nodes  arising 
from  an  artificial  element,  say  a  e  A,  into  a  source  node  s  and  a  destination 
node  d.  Let  the  arcs  out  of  s  in  Sp(I)  be  those  out  of  a  in  G(I)  and  the  arcs 
into  d  in  Sp(l)  be  those  into  a  in  G(I).  The  following  result  provides  the  basis 
for  an  algorithm  to  solve  (Pk). 

Theorem  2.  [BCG,  Theorem  5].  Let  I  be  optimal  for  (Pk|  H),  where  Echcf, 

and  let  Sp(i)  be  obtained  from  G(I)  by  splitting  an  artificial  element  a  into  s 
and  d. 

(i)  Problem  (Pk|  H  -  a)  has  no  solution  if  and  only  if  there  is  no  s-d 
dipath  in  Sp(I). 

(ii)  Let  P  be  a  shortest  s-d  dipath  such  that  every  s-d  dipath  defined 
on  a  subset  of  its  nodes  has  a  strictly  larger  weight,  and  let  £  and  £‘  be  the 
nodes  of  P  in  V  and  V,  respectively,  (s  and  d  both  give  rise  to  a  in  £.)  Then 
r  - 1  -  £  +  £'  is  optimal  for  (Pk|  H  -  a). 

We  use  this  to  solve  (Pk)  as  follows.  Start  with  H  -  F  and  the  set  I  -  A. 
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Construct  G(I)  and  then  Sp(l)  by  splitting  one  node  arising  from  an  artificial 
element  a.  Find  a  shortest  dipath  from  the  resulting  source  s  to  destination 
d,  with  the  added  condition  that  every  s-d  dipath  on  a  proper  subset  of  its 
nodes  has  a  strictly  larger  weight.  Then  the  set  I'  defined  in  Theorem  2(ii)  is 
optimal  for  (Pk|  H  -  a).  Repeat  this  process  with  r  in  place  of  I  and  H  -  a  in 

place  of  H,  until  all  artificial  elements  have  been  split.  When  this  occurs, 
we  have  an  optimal  solution  for  (Pk). 

We  close  this  section  with  a  result  that  will  be  useful  in  the  proof  of 
Theorem  6. 


Theorem  3.  [BCG,  Corollary].  Let  e  6  V  and  e'  <=  V\ 

(a)  If  P  is  an  e-e'  dipath  in  Sp(l)  and  I  - 1  +  r  is  a  dependent  set  in  M1 

then  there  is  an  e-e'  dipath  PQ  in  Sp(I)  such  that  the  intermediate  nodes  of  PQ 
constitute  a  proper  subset  of  those  of  P,  and  w(PQ)  £  w(P). 

(b)  If  P  is  an  e'-e  dipath  in  Sp(I)  and  I  -  2  +  r  is  a  dependent  set  in  M2 
then  there  is  an  e'-e  dipath  PQ  in  Sp(I)  such  that  the  intermediate  nodes  of  P0 
constitute  a  proper  subset  of  those  of  P,  and  w(P0)  £  w(P). 


3.  The  Algorithm 

In  this  section  we  show  how  to  solve  (P).  The  question  of  finding  a 
first  basis  B  is  not  addressed  until  later  in  the  section. 

For  convenience,  we  index  each  element  of  E  to  identify  the  subset  or 
state  of  the  partition  E1 ,  E2, ....  Em  to  which  it  belongs.  Thus,  e-,  and  fj  are 

elements  of  Ej  and  Ej,  respectively.  Also,  the  symbol '  is  used  to  denote 

whether  or  not  an  element  is  in  B,  so  that  ej  e  B  but  ej'  e  E  -  B. 

We  define  the  state  graph  of  B,  denoted  S(B)  -  (N,  A),  as  the 
following  directed  graph.  The  node  set  is  given  by  N  -  {v^ ,  v2, ....  vm},  where 

node  Vj  corresponds  to  state  i,  i  -  1 , 2 . m.  For  each  ej  in  B  we  construct  an 

arc  a(ej(  6j')  with  weight  w(a(6j,  ej'))  -  w(ej')  -  w(ej)  directed  from  node  Vj 

to  node  Vj  provided  (ej,  ©j*)  is  a  B-swap.  We  call  these  the  forward  arcs  of 

S(B).  Note  that,  in  general,  S(B)  is  a  multigraph  and  can  have  loops.  For  each 
6j'  in  E  -  B,  construct  a(6j',  6|)  from  Vj  to  Vj,  with  weight  0,  provided 

B  -  ej  +  6j'  does  not  violate  the  cardinality  conditions  of  E(  and  Ej.  We  call 

these  the  backward  arcs  of  S(B).  Note  that  the  existence  of  a  backward  arc 
from  Vj  to  Vj  indicates  that  we  can  exchange  any  element  e(  in  B  n  Ej  for  any 

element  6j'  in  (E  -  B)  n  Ej  without  violating  the  cardinality  conditions  on 

B  n  Ej  and  B  n  Ej.  All  such  backward  arcs  a(6j\  e()  appear  in  parallel  in  S(B). 

We  could  have  constructed  S(B)  using  G(B),  the  bipartite  auxiliary 
digraph  discussed  in  the  previous  section,  as  follows.  Let  M1  »  M  and  let  M2 

be  the  matroid  whose  bases  are  elements  of  the  set 

B*  -  {X  c  E:  lj  <;  |X  n  Ej|  <;  ujf  i  -  1 , 2, ....  m,  |X|  -  r}. 

Note  that  (Pr)  defined  by  M1  and  M2  is  equivalent  to  problem  (P).  Construct 
G(B).  Consider,  for  i  -  1 , 2, ....  m,  the  set  of  nodes  corresponding  to  Ej  and 
identify  this  set  as  the  node  v,.  Thus,  each  arc  (ej,  ej’)  or  (fj’,  fj)  becomes  an 
arc  joining  Vj  to  Vj,  and  this  new  graph  is  clearly  S(B). 

Let  C  be  a  dicycle  of  S(B).  In  this  paper,  dicycles  can  have  repeated 
nodes,  but  no  repeated  arcs.  Note  that  this  does  not  exclude  the  possibility 
of  repeated  elements  of  E.  We  say  that  [3  c  B  and  (3'  c  E  -  B  give  rise  to  C  if 


p  -  {6j  e  B  :  a(e,,  ep  is  a  forward  arc  of  C},  and 

p*  -  {©j*  e  E  -  B  :  a(ejf  ej')  is  a  forward  arc  of  C}. 

We  define  the  weight  of  this  dicycle,  w(C),  to  be  the  sum  of  the  weights  of 
the  arcs  in  the  dicycie.  Then  w(C)  -  w(P')  -  w(p),  where  p  c  B  and  p'cE-B 
are  the  sets  of  elements  which  give  rise  to  the  dicycle.  Dipaths  of  S(B)  have 
similar  definitions,  with  the  following  exception.  If  the  last  arc  of  a  dipath 
P  is  a  backward  arc,  define 

p  -  {ej  €  B  :  a(ej,  ep  is  a  forward  arc  of  P 

or  a(ek\  ep  is  the  last  arc  of  P}. 

Note  that,  if  p  c  B  and  p'cE-B  give  rise  to  an  element  disjoint 
dicycle  C  in  S(B),  then  they  also  give  rise  to  a  dicycle  of  the  form 
{a(ej,  ep,  a(ep  ek), ....  a(e,,  ep'),  a(eq\  ej)}.  (*) 

There  are  four  ways  in  which  C  could  fail  to  have  this  form:  C  could  contain, 
as  subsequences, 

(i)  consecutive  forward  arcs  {a(ej,  ep,  a(fj,  ek')}, 

(ii)  consecutive  backward  arcs  {a(ej\  ej),  a(fj',  ek)}, 

(iii)  a  forward-backward  pair  with  different  elements  from 

E  -  B,  {a(e|,  ep,  a(fp  ek)},  or 

(iv)  a  backward-forward  pair  with  different  elements  from  B, 
{a(ep  fk),  a(ek,  e,')}. 

In  case  (i)  the  loop  a(ej\  fj)  must  be  a  backward  arc  in  S(B)  since  swapping 
these  two  elements  in  B  will  not  change  |B  n  Ej|.  So  replace  the  consecutive 
forward  arcs  in  C  with  {a(ejt  ep,  a(ep  fj),  a(fj,  ek')}.  (Recall  that  we  only 

require  dicycles  to  be  arc-disjoint;  so  loops  are  permitted.) 

The  result  of  the  consecutive  backward  arcs  in  case  (ii)  is  that  the 
cardinality  of  Ej  increases  by  one  in  B  -  {ej,  ek}  +  {ej',  fj'},  with  the 

cardinality  of  Ek  decreasing  by  one.  Ej  realizes  no  change  in  cardinality  in 

from  this  string.  Thus,  a(ej\  ek)  must  be  a  backward  arc  in  S(B),  and  we  can 

replace  the  consecutive  backward  arcs  in  case  (ii)  with  the  single  backward 
arc  a(ej’,  ek). 

In  case  (iii),  the  backward  arc  a(fp  ek)  can  be  replaced  with  the 
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parallel  arc  a(6j\  ek).  Similarly,  a(ej\  fk)  can  be  replaced  with  a(ej\  ek)  in 
case  (iv). 

Since  these  steps  (which  result  in  a  dicycle  that  has  the  form  (*))  only 
affect  the  backward  arcs  in  the  dicycle,  p  and  p'  are  unchanged. 

The  following  result  is  immediate. 

Lemma  1.  There  is  a  one-to-one  correspondence  between  element-disjoint 
dicycles  in  S(B)  having  the  form  (♦)  and  node-disjoint  dicycles  in  G(B). 

Theorem  4.  Assume  B  is  feasible  for  (P). 

(i)  B  is  optimal  for  (P)  if  and  only  if  there  are  no  negative  dicycles  in 

S(B). 

(ii)  Let  C  be  a  negative  dicycle  in  S(B)  with  no  negative  dicycle  on  a 
subset  of  the  elements  p  +  p'  which  give  rise  to  C.  Then  B'  =  B  -  p  +  P'  is  a 
base  such  that  w(B')  <  w(B). 

Proof.  This  result  follows  directly  from  Lemma  1 ,  Theorem  1 ,  and  the 
following  observation.  If  a  negative  dicycle  in  S(B)  has  a  repeated  element, 
that  dicycle  induces  two  new  dicycles,  one  of  which  must  be  negative.  Thus, 
if  there  is  a  negative  dicycle,  there  is  a  negative  dicycle  without  repeated 
elements.  // 

Theorem  4  shows  the  equivalence  of  solving  (P)  and  finding  negative 
weight  dicycles  in  S(B)  with  the  property  that  no  negative  dicycle  exists  on  a 
subset  of  the  elements  which  give  rise  to  it.  This  gives  the  basis  for  a 
primal  algorithm:  start  with  a  feasible  B  and  use  S(B)  to  find  an  improved 
solution.  Continue  in  this  manner  until  we  find  an  S(B)  with  no  negative 
dicycles. 

We  briefly  discuss  the  complexity  issues  of  this  algorithm:  Firstly, 
one  must  be  able  to  find  a  negative  dicycle  in  S(B)  with  the  required 
property.  Also,  the  number  of  iterations  needed  to  reach  the  optimal  solution 
can  be  relatively  large:  even  though  each  iteration  yields  a  better  solution, 
the  improvement  may  be  slight.  For  these  reasons  the  complexity  is  high, 
and  we  do  not  pursue  our  investigation  of  this  primal  algorithm.  (We  do  note, 
however,  that  this  approach  may  have  merits  in  the  context  of  sensitivity 
analysis:  Given  an  optimal  solution  for  a  certain  set  of  weights, 


reoptimizing  on  a  perturbed  set  of  weights  using  this  method  may  be  more 
efficient  in  practice  than  starting  from  scratch,  as  other  algorithms  for 
matroid  intersection  would  require.) 

Instead  we  use  these  results  to  develop  a  dual  algorithm.  First,  using 
a  modification  of  the  greedy  algorithm,  we  can  find  an  initial  solution  BQ. 

Start  with  B0  -  0  and,  at  each  iteration,  add  an  element  ee  E  -  B0  to  B0  if 

(i)  Bq  +  e  is  independent  in  M, 

(ii)  e  has  the  smallest  weight  of  all  such  f  e  E  -  B0>  and 

(iii)  |(B0  +  e)  n  Ej|  £  Uj  for  i  -  1 ,  2, ...,  m. 

Stop  when  either  no  such  e  can  be  found,  |B0|  -  r,  or,  for  the  e  chosen, 

r  -  |B0  +  e|  <  X{max(lj  -  |(B0  +  e)  n  Ej|,  0):  i  -  1 , 2 . m}. 

Violating  this  last  condition  would  result  in  a  B0  which  could  not  satisfy 
both  |Bq|  -  r  and  all  lower  bound  conditions.  At  the  termination  of  this 
procedure  we  clearly  have  the  least  weight  independent  set  with  |B0| 
elements. 

If  care  has  been  taken  to  delete  loops  of  the  matroid  M  from  E,  B0  will 
have  at  least  one  element.  If  |B0|  » r,  (P)  is  solved.  Otherwise,  consider  a  set 
A  of  r  -  |B0|  artificial  elements,  where  E  n  A  -  0.  We  assign  to  each 
artificial  element  a  state  as  follows.  Extend  E1  by  unioning 
max(0,  l-|  -  |B0  n  E-jl)  artificial  elements.  Then  extend  E2,  E3 . Em 

similarly.  Any  remaining  artificial  elements  can  then  be  unioned  to  any  of 
the  states,  as  long  as  the  resulting  states  satisfy  |(B0  +  A)  n  Ej|  <  Uj.  Note 

that,  if  this  cannot  be  done,  (P)  is  infeasible. 

Let  F  -  E  +  A  and  define  a  matroid  M(F)  on  the  set  F  as  follows:  for 
I  c  E  and  J  c  A,  I  +  J  is  independent  in  M(F)  if  and  only  if  I  is  independent  in  M 
The  problem  (P)  relative  to  M(F)  will  be  denoted  (P|  F);  we  now  wish  to  find 
the  least  weight  independent  subset  B  of  F  such  that  I;  <  |B  n  Ej|  £  Uj,  for 

i  -  1 , 2, ....  m,  and  |B|  -  r.  Note  that,  if  we  define  the  matroids  M1  and  M2  as 
at  the  beginning  of  this  section  and  M-j(F)  and  M2(F)  in  the  obvious  way,  then 
(P|  F)  is  equivalent  to  (Pr|  F)  defined  by  M^F)  and  M2(F). 
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To  define  (P|  F)  completely  we  need  to  assign  weights  to  the  artificial 
elements.  If  we  give  each  artificial  a  large  positive  weight,  we  can 
guarantee  that,  if  (P)  is  feasible,  then  an  optimal  solution  to  (P|  F)  contains 
no  artificial  elements  and  is,  hence,  optimal  for  (P)  as  well.  Thus,  we  can 
use  Bq  +  A  as  a  starting  solution,  and  solve  (P|  F)  using  the  primal  algorithm 

whose  basis  is  given  earlier  in  this  section. 

However,  we  have  a  different  algorithm  in  mind.  We  will  give  large 
negative  weights  to  all  artificial  elements.  Thus,  the  initial  solution 
B  -  B0  +  A  will  be  optimal  for  (P|  F). 

Define  the  digraph  S*(B)  from  S(B)  by  creating  two  nodes,  a  source 
node  s  and  a  destination  node  d,  from  an  artificial  element,  say  a  e  A.  We 
also  allow  s  and  d  to  denote  the  resulting  elements.  Assume  a  was  assigned 
to  Ej.  Replace  arcs  of  the  form  a(a,  op  in  S(B)  with  arcs  a(s,  ej')  in  S*(B), 

and  replace  arcs  of  the  form  a(ej',  a)  in  S(B)  with  a(6|',  d)  in  S*(B). 

Note  that  we  could  have  constructed  S*(B)  by  forming  Sp(B)  from  G(B), 
and  then  shrinking  groups  of  nodes  (excluding  s  and  d)  corresponding  to  each 
of  the  states  into  the  nodes  v1 ,  v2,  ...,vm.  Thus  the  following  theorem  is  a 

direct  result  of  Theorem  2. 

Theorem  5.  Let  B  be  optimal  for  (P|  H),  where  EchcF,  and  let  S*(B)  be 
obtained  from  S(B)  by  creating  s  and  d  from  artificial  element  a. 

(i)  Problem  (P|  H  -  a)  has  no  solution  if  and  only  if  there  is  no  s-d 
dipath  in  S*(B). 

(ii)  Let  P  be  a  shortest  s-d  dipath  such  that  every  s-d  dipath  defined 
on  a  subset  of  the  elements  which  give  rise  to  P  has  a  strictly  larger  weight, 
and  let  p  and  P'  be  the  elements  of  B  and  H  -  B,  respectively,  that  give  rise  to 
the  nodes  of  P.  (s  and  d  both  give  rise  to  a  in  p.)  Then  B'  -  B  -  p  +  p’  is 
optimal  for  (P|  H  -  a). 

Theorem  5  provides  the  basis  for  our  algorithm.  Start  with  H  -  F  and 
B  -  Bq  +  A,  which  is  optimal  for  (P|  F).  Construct  S(B)  and  then  S*(B)  by 

creating  s  and  d  from  artificial  element  a.  Find  a  shortest  s-d  dipath  P  such 
that  every  dipath  on  a  subset  of  the  elements  which  give  rise  to  P  has  a 
strictly  larger  weight.  Then  the  set  B'  defined  in  Theorem  5(ii)  is  optimal 


for  (P|  H  -  a).  Repeat  this  process  with  B'  in  place  of  B  and  H  -  a  in  place  of 
H,  until  the  resulting  B'  contains  no  artificial  elements.  At  this  step,  B'  is 
optimal  for  (P). 

Note  that  the  arc  weights  of  S(B0  +  A)  are  all  nonnegative.  To  see  this, 

first  consider  arcs  of  the  form  a(a,  ep,  where  a  is  artificial.  Since  w(a) 

has  been  chosen  small  enough,  w(a(a,  ep)  -  w(ep  -  w(a)  >  0.  Now  consider 

an  arc  a(©j,  ep,  where  ej  €  B0.  Since  this  arc  is  in  S(B0  +  A),  we  have 

B0  -  ©j  +  6j'  independent  in  M.  But  then  w(ej)  £  w(ep  by  requirements  (i)  and 

(ii)  of  our  modified  greedy  approach  in  the  construction  of  B0;  so 

w(a(ej,  ep)  £  C.  Finally,  when  e  is  not  in  B0  +  A,  w(a(e,  f))  -  0. 

All  arc  weights  nonnegative  in  S(B0  +  A)  implies  the  same  for 

S*(B0  +  A);  so  no  shortest  s-d  dipath  will  repeat  nodes  (or,  hence,  elements). 

Thus,  for  all  pairs  of  nodes,  we  only  need  to  consider  the  shortest  of  the  arcs 
joining  Vj  to  Vj  when  searching  for  the  shortest  s-d  dipath,  and  we  can  apply 

Dijkstra's  algorithm  to  find  a  shortest  s-d  dipath  P.  Furthermore,  the 
requirement  that  every  dipath  on  a  subset  of  the  elements  which  give  rise  to 
P  has  a  strictly  larger  weight  can  be  obtained  by  adding  a  small  positive  e  to 
all  the  arc  weights  of  S*(B0  +  A). 

We  now  show  that  nonnegative  weights  on  the  arcs  of  S(B)  can  be 
preserved  throughout  the  algorithm.  Define  a  variable  D(Vj)  associated  with 

every  node  of  the  digraph  S*(B)  and  a  reduced  weight 
w'(a(ej,  ej'))  -  w(a(ej,  ep)  +  D(Vj)  -  D(vp 

associated  with  each  forward  arc  of  S(B).  Similarly,  for  each  backward  arc 
define 

w'(a(ek\  e,))  -  D(vk)  -  D(V|). 

Note  that,  in  terms  of  the  reduced  weights,  the  length  of  an  s-d  dipath  P  is 
w'(P)  -  w(P)  +  D(s)  -  D(d)  since,  for  any  intermediate  node  vj,  the  variable 

D(Vj)  cancels  out  on  the  two  arcs  of  P  that  contain  Vj.  Since  D(s)  -  D(d)  is  a 

constant  that  does  not  depend  on  P,  it  is  equivalent  to  solve  the  shortest 
dipath  problem  in  S*(B)  with  the  reduced  weights  w*  instead  of  the  original 
weights  w.  We  provide  a  choice  for  the  variables  D(Vj)  that  guarantees 
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nonnegative  reduced  weights  from  one  iteration  to  another  in  the  next 
theorem. 

Theorem  6.  Let  B  be  optimal  for  (P|  H)  where  ECHCF,  and  assume  that  an 
s-d  dipath  exists  in  S*(B).  Set  D(V|)  to  be  the  length  of  a  shortest  s-Vj 

dipath  in  S*(B),  for  i  -  1 , 2 . m.  Let  B'  be  as  defined  in  Theorem  5(ii).  Then, 

the  reduced  weights  W(a(ej,  @p)  -  w(a(ej,  ep)  +  D(v;)  -  D(Vj)  and 

w'(a(ek\  oj)) «  D(vk)  -  D(vj)  are  nonnegative  for  every  forward  arc  a(ejp  ep 

and  every  backward  arc  a(ek\  e()  of  the  digraph  S(B')  for  problem  (P|  H  -  a). 

It  was  noted  above  that  every  dicycle  in  S(B)  can  be  assumed  to  have 
the  form  (*).  Using  the  same  argument,  we  can  assume  every  s-d  dipath  in 
S*(B)  to  have  the  form 

P  -  {a(s,  ep,  a(ej\  ep,  a(ej,  ek‘) . a(et\  d)}.  (*♦) 

Note  that,  since  the  only  arcs  affected  by  changing  an  s-d  dipath  to  one  of 
the  form  (**)  are  backward  arcs,  putting  an  s-d  dipath  into  this  form  does  not 
change  the  weight  of  the  dipath. 

If  an  s-d  dipath  is  of  the  form  (**),  the  elements  which  give  rise  to 
P  are  all  elements  represented  in  an  arc  of  P  since  these  are  exactly  the 
elements  which  are  represented  in  the  forward  arcs  of  P.  For  the  proof 
of  Theorem  6  we  assume  that  the  shortest  s-d  dipath  has  been  put  in  the 
form  (**).  We  also  need  the  following  remark,  which  is  a  direct  result  of 
Theorem  3. 

Remark  3.  Let  P  be  a  shortest  s-d  dipath  in  S*(B)  such  that  every  s-d 
dipath  of  the  form  (**)  on  a  subset  of  the  elements  which  give  rise  to  P  has 
strictly  larger  weight,  and  let  P0  be  some  subpath  of  P  from  s  to  another 

node.  Let  p0  and  0O’  be  the  sets  of  elements  which  give  rise  to  P0. 

(a)  If  the  last  arc  on  P0  is  a  forward  arc,  then  B  -  (30  +  jJ0'  is 
independent  in  M. 

(b)  If  the  last  arc  on  PQ  is  a  backward  arc,  then  B  -  p0  +  p0'  satisfies 
the  cardinality  conditions  on  each  Ej. 


Proof  of  Theorem  6.  First  we  show  that  the  choice  for  D(Vj)  gives 
nonnegative  reduced  weights  on  the  arcs  of  S(B).  Consider  a(ejf  ep  in  S(B), 
and  suppose  w(a(ej,  ej'))  +  D(Vj)  -  D(Vj)  <  0.  This  gives 
D(vj)  +  w(a(ei,ej,))<D(Vj), 

which  implies  that  the  shortest  s-Vj  dipath  together  with  the  arc  a(ejt  e  ’) 
gives  a  shorter  distance  from  s  to  Vj  than  D(Vj),  a  contradiction.  Similarly, 
for  each  backward  arc  a(fj\  fj)  in  S(B),  w(a(fp  fj))  +  D(Vj)  -  D(Vj)  £  0.  Thus, 
our  choice  of  D(vj)  gives  nonnegative  reduced  weights  on  the  arcs  of  S(B). 

Let  P  be  the  shortest  s-d  dipath  that  gives  rise  to  B'  in  Theorem  5(ii), 
and  let  Pj  be  the  subpath  of  P  from  s  to  some  node  Vj.  Let  p-j  (i)  and  p-j'fl)  be 

the  sets  of  elements  of  B  and  E  -  B,  respectively,  represented  by  the  forward 
arcs  of  Pj,  and  let  p2(i)  and  P^O)  be  the  corresponding  sets  of  elements 

represented  by  the  backward  arcs  of  Pj.  By  Remark  4,  B1  (i)  -  B  -  p1  (i)  + 

P-l'O)  is  independent  in  M,  and  B2(i)  -  B  -  p2(i)  -1-  p2'(i)  satisfies  the 
cardinality  conditions  on  each  Ej.  For  ease  of  notation  we  allow  v;  to 
represent  s  and  d,  as  well  as  the  nodes  representing  states  E1 ,  E2, ...,  Em. 

'.Ve  construct  a  digraph  Nj  as  follows.  The  node  set  of  Nj  is  that  of 
S(B).  If  B1  (i)  -  ej  +  ek'  is  independent  in  M,  then  a(ej,  ek')  is  a  forward  arc  of 
Nj;  if  B2(i)  -  ej  +  ek'  satisfies  the  cardinality  conditions  of  each  Ej,  then 
a(ek\  ej)  is  a  backward  arc  of  Nj.  Note  that,  in  the  latter  case,  we  call 
(ek‘,  op  a  B2(i)-swap  with  respect  to  the  partition  matroid. 

We  will  show  that  this  choice  of  D(Vj)  gives  nonnegative  reduced 
weights  on  the  arcs  of  Nj  by  induction  on  the  nodes  of  P,  starting  from  Vj  =  s. 
Note  that  when  we  reach  Vj  -  d,  we  will  have  Nj  -  S(B')  for  (P|  H  -  a),  and  the 
theorem  will  be  proved. 

When  Vj  -  s,  Nj  -  S(B),  and  the  result  has  been  proved  above.  Let  Vj  be 
node  of  P  for  which  the  result  holds,  and  let  Vj  be  the  node  following  Vj  in  P. 
There  are  two  cases  to  consider,  as  the  arc  joining  Vj  to  Vj  can  be  a  forward 
arc  or  a  backward  arc. 


Suppose  the  arc  joining  v;  to  Vj  in  P  is  a  forward  arc  a(6j,  ep.  In  this 

case  the  backward  arcs  of  Nj  are  the  same  as  those  of  Nj;  hence,  they  have 

nonnegative  reduced  weights  by  induction.  Thus,  we  only  need  to  show 
D(Vq)  -  D(Vp)  £  w(eq’)  -  w(ep)  for  every  B1  Q-swap  (ep,  6q').  Note  that  if 

(ep,  6q')  is  also  a  B1  (i)-swap  we  are  done.  So  only  consider  (j)-swaps 

which  are  not  B-j  (i)-swaps. 

If  (ep.  ®q’)  -  (®j',  ®j).  then 

D(Vj)  -  D(vj)  +  w(ep  -  w(6j)  (t) 

since  a(ej,  ep  is  on  the  shortest  s-d  dipath  P,  and  the  desired  inequality  is 
immediate. 

For  the  remaining  cases,  consider  B*  -  B-j  (i)  -  ej  +  ej'  -  ep  +  which 
is  independent  in  M.  Suppose  ep  -  ©j’  and  eq’  *  ej.  Then  B*  -  B-|  (i)  -  ej  +  6q' 
and  (ej,  e^)  is  a  Bj-swap.  Therefore,  w(eq')  -  w(ej)  +  D(Vj)  -  D(vq)  >0  by 
induction.  Combining  this  inequality  with  (f)  and  ep  -  ej',  we  get 
w(eq')  -  w(ep)  +  D(vp)  -  D(Vq)  £  0,  and  the  result  follows.  A  similar  argument 
applies  when  ep  *  ej'  and  -  ej. 

Finally,  suppose  ep  *  ej'  and  Oq  *  ej.  Since  (ep,  eq')  is  not  a 
B1  (i)-swap,  yet  B*  is  a  base,  it  follows  from  Remarks  1  and  2  that  both 
(ep,  ej')  and  (ej,  eq’)  are  B1  (i)-swaps.  Thus,  w(ep  -  w(ep)  +  D(vp)  -  D(vj)  £  0, 
and  w(eq')  -  w(ej)  +  D(ej)  -  D(eq)  £  0.  Summing  these  inequalities  and 
applying  (t),  we  get  w(eq’)  -  w(ep)  +  D(Vp)  -  D(vp)  ^  0. 

Now  suppose  the  arc  joining  Vj  to  Vj  in  P  is  a  backward  arc  a(fj\  fj). 
Then  B1 G)  -  B1  (i),  and  the  forward  arcs  of  Nj  are  the  same  as  those  of  Nj. 

Hence,  they  have  nonnegative  reduced  weights  by  induction.  Thus,  we  only 
need  to  show  D(vp)  -  D(vp)  £  0  for  every  B2(j)-swap  (ep\  ep).  Note  that  if 

(ep\  Oq)  is  also  a  B2(i)-swap  we  are  done  by  induction.  So  only  consider 

B2(j) -swaps  which  are  not  B2(i)-swaps. 

If  (®p.  ®q)  -  (fj,  fj'),  then 
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since  a(fj',  fj)  is  on  the  shortest  s-d  dipath  P,  and  the  desired  inequality  is 
immediate. 

For  the  remaining  cases  consider  B“  -  B2(i)  -  fj  +  fj'  -  ep  +  ep',  which 
satisfies  the  cardinality  conditions  for  each  Ej  since  (ep\  eq)  is  a 
B2(j)-swap.  Suppose  ep'  -  fj  and  eq  *  fj'.  Then  B**  ■  B2(i)  -  eq  +  fj',  and 
(fj',  eq)  is  a  B2(i)-swap.  Hence,  D(vj)  -  D(vq)  £  0  by  induction.  This  and 
equation  (ft)  show  that  D(Vp)  -  D(Vq).  The  proof  is  similar  when  ep'  *  fj  and 

eq  *  fi'- 

Rnally,  suppose  ep'  *  fj  and  eq  *  fj'.  Since  (ep\  eq)  is  not  a 
B2(i)-swap,  yet  B**  is  independent  in  the  partition  matroid,  it  follows  from 
Remarks  1  and  2  that  both  (ep'(  fj)  and  (fj',  ep)  are  B2(i)-swaps.  Thus, 
induction  yields  D(vp)  -  D(Vj)  >  0  and  D(Vj)  -  D(Vq)  >  0.  Summing  these 
inequalities  and  applying  (tt),  we  get  D(vp)  -  D(vq)  £  0.  // 

We  now  summarize  the  complexity  of  the  algorithm.  O(r)  iterations 
will  be  needed,  one  iteration  for  the  removal  of  each  artificial  element  in 
the  initial  solution.  In  each  iteration  we  need  to  construct  S(B),  find  a 
shortest  s-d  dipath,  and  update  the  variables  D(Vj).  Recall  that  |E|  -  n,  and  c 

represents  the  time  required  to  find  the  circuit  of  I  +  e’  in  M  (or  show  that 
none  exists),  where  I  c  E  is  independent. 

To  construct  S(B)  for  the  problem  (P|  H),  E  c  H  c  F,  we  solve  the 
following  circuit  recognition  problem  for  each  e'  e  H  -  B:  Check  whether 
B  +  e'  is  dependent  in  M(H),  and,  if  it  is,  find  the  unique  circuit  of  B  +  e'.  This 
task  requires  time  c  for  each  e’  e  H  -  B.  So  the  complexity  of  constructing 
S(B)  is  O(nc).  In  some  instances  it  is  possible  to  speed  up  the  construction 
of  S(B),  as  not  all  parallel  arcs  need  to  be  constructed,  but  only  the  shortest. 
An  example  of  this  will  be  given  in  Section  5. 

Using  the  fast  implementation  of  Fredman  and  Tarjan  (see  [FT]),  the 
complexity  of  finding  a  shortest  s-d  dipath  by  Dijkstra’s  algorithm  is  at  . 
most  0(p  +  m  log  m),  where  p  is  the  cardinality  of  the  set  of  arcs  from  S*(B) 
used.  Recall  that  we  do  not  need  loops  and  only  consider  the  shortest  in  each 
group  of  parallel  arcs  to  find  a  shortest  s-d  dipath. 


Note  that  the  variables  D(vj)  needed  in  Theorem  6  are  actually 

computed  in  the  course  of  finding  a  shortest  s-d  dipath.  Thus,  no  extra 
computations  are  needed.  Also,  nc  will  be  larger  than  p.  Therefore,  the 
overall  complexity  of  the  algorithm  is  0(r(nc  +  m  log  m)). 


4.  Weighted  bipartite  matching. 

Consider  a  graph  G  -  (V,  E)  and  a  weight  function  w:  E  ->  R.  The 
matching  problem  consists  of  finding  a  minimum  weight  subset  F  c  E  such 
that  each  node  of  V  is  incident  with  exactly  one  member  of  F.  In  this  section 
we  assume  that  G  is  a  bipartite  graph  where  the  partition  V1  u  V2  -  V  is 

such  that  |V1 1  -  |V2|.  The  weighted  bipartite  matching  problem,  also  called 

the  assignment  problem,  is  a  special  instance  of  problem  (P)  where  the 
partition  of  E  is  induced  by  the  nodes  of  V2.  Specifically,  e  e  Ej  if  and  only 

if  it  is  incident  with  node  v-,  c  V2.  For  each  i,  lj «  Uj » 1 .  The  matroid  M  is 

also  a  partition  matroid,  namely  the  one  induced  by  the  node  set  V1 . 

For  this  problem  r  -  |V|/2,  n  -  |E|,  and  m  -  |V|/2.  To  obtain  the 
complexity  of  our  algorithm  for  this  problem,  we  have  to  determine  the 
complexity  of  constructing  S(B).  Note  that  the  circuits  of  M  have  length  2; 
so  the  number  of  forward  arcs  of  S(B)  is  at  most  |E|.  Similarly,  there  are  at 
most  |E|  backward  arcs. 

Thus,  the  construction  of  S(B)  requires  time  0(|E|).  Since  the  Fredman 
Tarjan  algorithm  finds  the  minimum  length  dipath  in  time  0(|E|  +  |V|  log  |V|), 
the  overall  complexity  is  0{|V|  |E|  +  |Vp  log  |V|). 

Note  that  our  algorithm  is  closely  related  to  the  usual  augmenting  path 
algorithm.  Although  this  algorithm  and  the  associated  bound  are  not  new,  it 
is  interesting  that  they  are  obtained  using  the  state  graph  approach.  Recall 
that  the  general  matroid  intersection  algorithm  only  gives  the  bound 
0(|V|2|E|). 
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5.  Minimum  spanning  trees. 

Let  G  -  (V,  E)  be  an  undirected  graph,  and  let  each  edge  of  G  be  assigned 
a  weight.  Let  V*  -  u2, ....  ^m.\)  be  a  stable  set  of  G;  that  is,  no  edge  of  E 

has  both  endpoints  in  V*.  For  i  -  1 , 2, ....  m-1 ,  let  Ef  denote  the  edges  which 

are  incident  with  node  \)j  in  G,  and  let  Em  be  the  remaining  edges.  For  each  i 

assign  two  integers  lj  and  Uj,  where  lj  £  Uj.  The  minimum  spanning  tree 

problem  with  restrictions  on  the  number  of  edges  in  each  Ej  is  an  instance  of 

(P):  M  is  the  graphic  matroid  defined  on  the  edge  set  E,  and  the  partition 
matroid  is  given  by  the  restriction  that  the  number  of  edges  from  state  i 
belongs  to  the  interval  [lj,  Uj]. 

An  interesting  feature  of  this  instance  is  that,  by  using  the  special 
structure  of  the  tree,  we  do  not  need  to  construct  the  entire  state  graph  S(B). 
However,  our  method  requires  that  we  work  with  spanning  trees.  So,  instead 
of  simply  extending  the  graphic  matroid  M  to  M(F)  as  described  in  Section  3, 
create  an  artificial  edge  in  G  for  each  artificial  element  so  that  the 
resulting  first  solution,  B0  +  A,  is  a  spanning  tree  in  the  extended  graph 

G*  -  (V,  E  +  A). 

Define  the  star  of  a  node  to  be  the  set  of  edges  incident  with  that 
node.  Suppose  the  nodes  are  ordered  from  1  to  |V|.  For  node  1  define  its  net 
star  to  be  its  star.  Then,  for  i  -  2, 3, ....  |V|,  define  the  net  star  of  node  i  to 
be  its  star,  excluding  edges  which  appear  in  the  stars  of  nodes  1,2 . i-1 . 

Let  V*  be  nodes  1  to  m-1 ,  and  let  the  remaining  nodes  be  nodes  m 
to  |V|.  Then  the  net  star  for  each  node  \>j  in  V*  is  Ej,  and  Em  consists  of  the 

net  stars  for  V  -  V*.  Begin  by  sorting  the  edges  in  each  net  star  in  ascending 
order  by  weight.  At  worst,  this  work  is  0(|E|  log  |V|).  Note  that  we  need  to 
sort  the  edge  weights  only  once  since  the  reduced  weights  only  add  a 
constant  to  all  the  edges  in  the  same  net  star.  At  each  iteration  of  our 
algorithm,  in  the  construction  of  each  new  state  graph,  proceed  as  follows: 

1 .  Set  up  (or  update)  the  current  basis  tree  with  a  root  and  predecessor 
indexing. 


V 
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2.  For  i  -  1 , 2 . |V|,  start  with  all  nodes  unlabeled,  and  perform  the 

following  labeling  procedure,  examining  nonbasic  edges  in  node  i's  net  star  in 
ascending  weight  order. 

(a)  Label  all  nodes  in  the  predecessor  path  from  node  i  to  the  root 
with  an  x,  and  label  i  with  a  y. 

(b)  Let  j  be  the  other  endpoint  of  the  nonbasic  edge  (i,  j)  currently 
under  examination.  Trace  the  predecessor  path  for  j  until  reaching  a 
labeled  node  k,  where,  possibly,  k  -  j.  The  edges  between  k  and  j  in  the 
tree  have  never  been  examined  for  the  current  node  i.  Check  each  as  an  edge 
to  swap  with  (i,  j),  either  to  create  a  forward  arc  in  the  state  graph,  or  to 
replace  an  existing  parallel  one  of  larger  weight.  As  this  check  occurs,  label 
the  nodes  between  nodes  k  and  j  in  the  tree,  excluding  k,  with  a  y. 

(c)  If  node  k  has  a  y  label,  repeat  (b)  with  the  next  edge  in  node  i's 
net  star.  Otherwise,  k  has  an  x  label  and,  hence,  lies  on  the  predecessor 
path  from  node  i  to  the  root.  In  this  case,  trace  from  node  k  to  node  i, 
and  consider  swapping  each  edge  encountered  with  (i,  j),  stopping  when 
encountering  a  node  with  a  y  label.  When  this  occurs,  label  each  node 
encountered  with  a  y. 

Note  that  all  nodes  of  every  cycle  formed  with  an  added  edge  will  be 
labeled  with  a  y,  although  possibly  not  all  of  the  cycle  will  be  examined  on 
the  current  step,  since  some  edges  may  have  been  examined  before.  Also,  at 
every  step,  the  edges  examined  will  be  exactly  the  edges  of  the  current  cycle 
that  have  never  been  examined  before.  Finally,  because  we  examine  the 
nonbasic  edges  of  each  net  star  in  order  by  increasing  weight,  we  find  the 
minimum  weight  arc  in  S(B)  from  vf  to  Vj  for  each  pair  (Vj,  Vj)  of  states. 

We  now  consider  the  complexity  of  our  algorithm.  During  the  pass 
through  the  net  star  for  node  i,  which  has  at  most  |V|  -  i  edges,  the  work  of 
examining  the  edges  to  be  dropped  will  be  on  the  order  of  the  number  of 
nodes  in  the  tree,  not  previously  labeled,  which  are  now  labeled.  So,  0(|V|) 
work  is  required  for  each  nonempty  net  star,  which  yields  0(| V|2)  work  total. 

Note  that  the  complexity  due  to  the  Fredman-Tarjan  shortest  dipath 
algorithm,  at  most  0(m2),  is  dominated  by  the  0(|V|2)  required  to  construct 
each  state  graph.  The  graphic  matroid  has  rank  at  most  |V|  - 1 .  Hence,  the 
overall  complexity  is  OflVr). 


In  the  special  case  where  lj  -  uf  -  2  for  i  -  1 , 2 . m-1 ,  and  lm  *  0  and 

um  -  |V|  -  2m  +  1 ,  the  greedy  phase  of  finding  a  starting  solution  B0  can 

substantially  reduce  the  number  of  iterations  needed.  For  i  -  1 , 2 . m-1 , 

add  a  constant  c,-  to  the  weights  of  the  edges  incident  with  node  uf  in  V*. 

We  assign  large  enough  values  to  the  constants  so  that  the  greedy  phase  will 
first  choose  |V|  -  2m  + 1  edges  in  Em  to  construct  the  spanning  tree.  Let  w 

be  the  weight  of  the  last  edge  added.  For  i  -  1 , 2, ....  m-1 ,  we  adjust  the 
constant  cf  so  that  the  minimum  weight  edge  incident  with  has  weight  w 

and,  therefore,  can  be  chosen  next  in  the  greedy  phase. 

This  extended  greedy  solution  B0  contains  |V|  -  m  edges,  so  only  m-1 

artificial  elements  are  needed.  They  can  be  chosen  so  that  B0  +  A  is  a  tree 

having  degree  2  at  each  of  the  nodes  of  V*.  Furthermore,  if  the  weights  of 
the  artificial  elements  are  chosen  small  enough,  B0  +  A  is  an  optimum 

solution  of  (P|  E  +  A)  and,  therefore,  a  valid  start  for  our  dual  algorithm. 

Hence,  only  m  - 1  iterations  of  the  algorithm  are  needed.  This  yields  an 
algorithm  of  complexity  0(m  |V|2). 

In  this  case  (P)  is  the  problem  of  finding  the  minimum  weight  spanning 
tree  in  G  such  that  every  node  in  V*  is  adjacent  to  exactly  two  edges  of  the 
tree.  Note  that  such  a  spanning  tree  is  a  relaxation  of  a  hamiltonian  path, 
namely,  a  simple  path  which  contains  each  node. 

Consider  the  following  variant  of  spanning  trees.  A  1-tree  is  a  graph 
with  nodes  {1 , 2, ...,  |V|}  consisting  of  a  spanning  tree  on  {2,  3, ....  |V|}  together 
with  two  edges  incident  with  node  1 .  The  traveling  salesman  problem  on 
G  -  (V,  E)  seeks  a  minimum  weight  tour,  a  cycle  which  passes  through  each 
node  exactly  once.  Observing  that  a  tour  is  precisely  a  1-tree  in  which  each 
node  has  degree  2,  Held  and  Karp  explored  approaches  to  the  traveling 
salesman  problem  which  involve  1-tree  relaxations  (see  [HK1]  and  [HK2]). 

Our  technique  permits  a  generalization  of  Held  and  Karp's  method  to 
yield  a  constrained  1-tree  as  follows.  Suppose  the  nodes  are  indexed  so  that 

{1 , 2, ....  m-1}  is  a  stable  set  of  G,  and  let  V*  -  {2,  3 . m-1}.  Then,  instead 

of  using  a  general  spanning  tree  on  {2,  3, ...,  |V|},  find  a  minimum  weight 
spanning  tree  such  that  every  node  in  V*  has  degree  2.  Add  the  two  edges 
incident  with  node  1  of  minimum  weight.  Clearly,  this  constrained  1-tree  is 
still  a  relaxation  of  a  minimum  weight  tour.  It  is  tighter  than  the  1-tree 


relaxation,  which  is  obtained  when  V*  -  0. 

A  still  tighter  relaxation  can  be  obtained  using  the  net  star 
construction.  For  each  node  i,  let  r;  be  the  number  of  edges  deleted  from  its 

star  to  create  its  net  star.  Then  we  can  stipulate  that  the  spanning  tree  on 
(2,  3 . |V|}  must  contain  at  least  max{0,  2  -  r{}  and  at  most  2  edges  from 

the  net  star  of  each  node  i  >  2,  which  introduces  additional  constraints  to  the 
requirement  that  every  node  in  V*  has  degree  2.  This  is  a  valid  relaxation  of 
the  traveling  salesman  problem  even  if  {1 , 2, ....  m-1}  is  not  a  stable  set. 

Note  that  the  complexity  of  the  algorithm  in  this  case  is  still  0(m  |V|2). 

Finally,  consider  problem  (P)  where  M  is  a  graphic  matroid  and  E1 ,  E2> 
....  Em  is  a  general  partition  of  the  edge  set.  Net  stars  can  still  be  used 
within  each  Ej;  therefore,  the  least  weight  arc  from  state  i  to  each  of  the 

other  states  can  be  found  in  0(|V|2).  This  yields  a  complexity  of  0(m|V|2)  for 
the  construction  of  the  state  graph  and  0(m|V|3)  for  the  overall  algorithm. 

This  improves  on  the  complexity  of  the  general  matroid  intersection 
algorithm  when  |E|  >  0(m|V|).  We  leave  open  the  question  of  improving  this 
complexity  bound. 
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